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INTRODUCTION 


The  present  preliminary  report,  by  K.  0.  Priedrichs, 
was  prepared  under  Navy  Contract  N6ori-201,  Task  Order 
No.  1.   It  contains  an  approximate  treatment  of  inter- 
actions of  continuous  motion  in  a  fluid  with  shock  waves 
of  moderate  strength.   In  particular,  the  general  method 
developed  is  applied  to  the  problem  of  formation  of  shocks 
in  a  compression  wave,  to  the  decay  of  shocks  under  the 
influence  of  rarefaction  waves,  and  to  the  propagation  of 
N -waves  and  other  flow  patterns  as  they  occur  in  supersonic 
flow  past  airfoils. 

Apart  from  its  theoretical  interest  the  method  leads 
to  simple  procedures  for  numerical  calculation. 

In  the  appendices  of  a  more  detailed  publication  the 
mathematical  theory  will  be  amplified,  further  applications, 
especially  to  the  theory  of  water  waves  in  shallow  channels 
will  be  presented,  and  the  accuracy  of  the  approximate 
results  will  be  demonstrated. 

R .  C  our  ant 

Director,  Institute  for 

Mathematics  and  Mechanics 
New  York  University 


FORMATION  AND  DECAY  OP  SHOCK  WAVES 
By 
K.  0.  Friedrichs 


To  determine  the  interaction  of  shocks  with  continuous 
gas  flow  is  a  difficult  problem  which,  however,  can  be 
handled,  with  great  labor,  by  the  method  of  finite  differ- 
ences.  If  the  shock  is  weak  or  of  medium  strength  the 
treatment  can  be  greatly  simplified  in  good  approximation, 
as  Chandrasekhar  has  observed.  "'   He  obtained  an  explicit 
approximate  solution  of  the  interaction  problem  in  a  special 
case.   It  will  be  shown  in  the  present  paper  that  a  slightly 
modified  approach  enables  one  to  obtain,  quite  generally, 
explicit  analytic  solutions  of  such  interaction  problems  in 
good  approximation.   This  approach  opens  up  the  possibility 
of  handling  analytically  a  number  of  interaction  problems 
for  one -dimensional  unsteady  and  two-dimensional  steady  gas 
flow,  which  so  far  have  been  inaccessible. 

1.   One  such  problem  is  that  of  the  formation  of  a  shock 
wave.  Suppose  a  cylindrical  tube  full  of  gas  is  closed  at  one 
end  by  a  piston  which  is  pushed  into  the  gas  with  finite  ac- 
celeration. We  know  that  a  simple  condensation  wave  is 
propagated  into  the  gas  and  that  at  some  later  time  this  wave 
leads  to  a  discontinuity  which  is  resolved  by  the  development 
of  a  shock.   The  strength  of  this  shock  begins  with  the  value 
zero  and  builds  up  gradually.   Therefore  any  method  good  for 

^"'S.  Chandrasekhar,  On  the  decay  of  plane  shock  waves. 
Aberdeen  Proving  Ground,  Ballistics  Research  Laboratory, 
Report  No.  423,  8  November  1943. 
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treating  weak  or  medium  shocks  should  be  applicable  at  early 
and  intermediate  stages  of  the  process.  Let  us  describe  the 
motion  in  an  (x,t) -diagram 
and  assume  the  acceleration 
to  be  constant  up  to  a 
certain  time  and  then  zero 
again.  The  piston  is  then 
represented  by  a  section  of 


Formation  of  shock  in  wave 
created  by  piston  moving 
with  constant  acceleration. 


1   2 
a  parabola  x  =  -g  bt  ,  con- 
tinued as  a  tangent  from 
some  point  on.   The  result- 
ing flow,  a  "simple"  com- 
pression wave,  may  be  de- 
scribed with  the  aid  of  the 
characteristic  lines  issuing 
from  the  piston  curve.  They 

are  straight  lines  and  carry  constant  values  of  all  pertinent 
quantities.   These  straight  lines,  continued  indefinitely, 
form  an  envelope  with  a  cusp  at  some  point  on  the  first 
characteristic  x  =  c  t  issuing  from  the  point  x  =  0,  t  =  0 
(cf.  Manualv"  ,  III  25,  p.  60).  The  cusp  region  between  the 
two  branches  of  the  envelope,  of  which  one  is  a  section  of 
the  first  characteristic  x  =  cQt,  is  covered  three  times  by 
the  straight  lines.  The  continuation  of  the  flow  from  the 
left  or  from  the  right  of  the  cusp  would  lead  to  a  different 
description.  This  contradiction  is  resolved  by  the  appearance 
of  a  shock  across  which  the  pertinent  quantities  of  the  flow 
undergo  a  jump.   The  shock  influences  the  flow  behind  it; 
more  specifically,  the  flow  differs  from  the  simple  compression 
wave  described  in  an  increasing  zone  represented  by  the  region 
in  the  (x,t) -plane  obtained  by  drawing  the  "cross-character- 
istic" backward  from  the  cusp  to  the  piston  curve •   (This  in- 


WL 


Supersonic  Plow  and  Shock  Waves,  For  the  Applied  Mathema- 
tics Panel,  National  Research  Defense  Committee,  by  the 
Applied  Mathematics  Group,  New  :'ork  University,  1944, 
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fluence  of  the  shock  on  the  flow  can  to  a  certain  extent 
be  described  in  terms  of  a  "reflected  wave",)   The  forma- 
tion of  this  shock  can  now  easily  be  described  on  the 
basis  of  the  method  to  be  presented.  "'' 

2.  The  problem  of  the  decay  of  a  shock  wave  is  en- 
countered when  the  piston,  through  whose  motion  a  shock 
was  formed  in  the  tube,  is  arrested  or  retracted.  The  rare- 
faction wave  sent  out  during  the  process  of  deceleration 
follows  the  shock  front  with  sound  speed;   hence  it  even- 
tually catches  up  with  it,  and  diminishes  its  strength. 
Let  us  assume  that  the  piston  was  instantaneously  acceler- 


ated to  constant  velocity; 
then  a  constant  shock  re- 
sults. Arresting  the 
piston  produces  a  simple 
rarefaction  wave.   (The 
total  phenomenon:   shock 
front  followed  by  a  rare- 
faction wave,  which  is 
modified  by  the  shock,  is 
often  called  "shock  wave".) 
The  width  of  an  undisturbed 
simple  rarefaction  wave  at 
any  time  increases  linearly 


Shock  wcjt-e-  Zone  c#Z) 


"X 


Pig.  2 


Decay  of  shock  wave  due 
to  arrest  of  piston. 


with  time.   The  shock  front,  through  its  diminishing  velocity, 
reduces  the  rate  of  increase  of  the  width  of  the  wave  zone. 
Whether  or  not  the  width  of  this  "decaying"  shock  wave  zone 
actually  decreasec  and  whether  or  not  the  interaction  process 
comes  to  an  end  at  a  finite  time  has  often  been  asked.  The 
approximate  treatment  to  be  presented  gives  a  definite 
(though  not  rigorously  proved)  answer:   the  width  of  the 


TO 


v.  Neumann  indicated  in  one  of  his  reports  some  years 
ago  that  Calkin  made  a  study  of  the  formation  problem, 
but  no  report  about  this  study  is  available. 
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decaying  shock  wave  zone  increases  like  the  square  root, 
and  the  strength  of  the  shock  decreases  like  the  reciprocal 
of  the  square  root,  of  the  time  elapsed  since  the  start  of 
the  process. 

Other  problems  of  formation  and  decay  of  shocks,  such. 
as  the  problem  of  "N-waves"  (Sec.  7)  and  the  problem  of  the 
decay  of  a  two-dimensional  steady  shock  (Sec.  10),  will  be 
described  later  on. 

3.  At  first  we  shall  present  the  method  of  treatment 
to  be  employed.  The  method  is  based  on  the  fact  that  the 
shock  transition  relation  agrees  in  second  order  with  the 
simple  wave  transition  relations.  More  specifically:   con- 
sider forward  facing  shocks  of  varying  strength  with  the 
same  state  ahead  of  them  and  expand  all  pertinent  quantities 
behind  the  wave  in  powers  of  any  quantity  that  can  serve  as 
a  measure  of  shock  strength,  such  as  the  velocity  difference* 
Do  the  same  with  the  quantities  behind  a  forward  facing 
simple  wave  having  the  same  state  ahead  of  it  as  the  shock. 
The  expansions  of  sound  speed  c  and  entropy  s  with  respect 
to  powers  of  the  strength  then  agree  for  the  shock  and  the 
simple  wave  up  to  terms  of  second  order;  the  same  is  true 
for  the  expansion  of  pressure  and  density.  Since  the  entropy 
is  constant  across  a  simple  wave,  this  statement  implies  the 
well  known  fact  that  the  change  of  entropy  across  a  shock  is 
of  third  order.  Since  the  Riemann  invariant  (see  below),  is 
constant  across  a  simple  wave,  the  change  of  the  Riemann 
invariant  across  a  shock  is  also  of  third  order.  This  state- 
ment deserves  considerable  emphasis,  which  it  has  not  received 
in  the  literature  where  it  was  occasionally,  but  somewhat 
casually,  mentioned.   (For  its  derivation  see  Appendix  BI.) 

The  method  to  be  presented  consists  of  two  steps.   In 
the  first  step  we  replace  two  shock  transition  relations  by 
the  corresponding  transition  relations  for  simple  waves. 
Specifically,  we  require  that  the  entropy  and  the  appropriate 
Riemann  invariant  do  not  change  across  the  shock.  Owing  to 
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this  simplification  there  is  no  back  reaction  of  the  shock 
on  the  flow  behind  it,  as  we  shall  justify  in  detail  below. 
Consequently,  we  may  describe  the  flow  in  the  following 
manner.  We  first  determine  the  simple  wave  independently 
of  the  shock  and  then  draw  the  shock  line  in  the  wave  region 
in  an  appropriate  manner. 

Suppose  at  the  time  t  =  0,  the  velocity  u  and  sound 
speed  c  are  given  as  functions  of  x  so  as  to  satisfy  the  re- 
lations 

(1)  ^u(x)  -  (1  -  JjF)    c(x)  -  const. 

for  a  forward  facing  simple  wave  (cf.  Manual  111-20).  Here 

2    Y— 1 

and  in  the  following  we  use  the  abbreviation  u     =  ^ty  • 

The  left  member  of  equation  (1)  is  just  the  Riemann  invariant 
mentioned  above.  At  the  point  x  =  0  the  quantities  u  and  c 
shall  suffer  a  discontinuity.  We  characterize  by  the  sub- 
scripts (o)  and  (1)  the  states  at  the  positive  and  negative 
side  of  the  cross  section  x  =  0.  The  requirement  that  the 
discontinuity  involves  a  compression  and  that  the  gas  crosses 
this  cross  section  from  the  positive  to  the  negative  side  is 
then  expressed  through  the  conditions 


U^  <  Ut   ,    C   <  C-,   • 


The  discontinuity  is  therefore  a  forward  facing  shock;  the 
shock  relations  are  satisfied  up  to  the  second  order  in  the 
shock  strength. 

We  now  determine  a  simple  wave  flow  for  t  >  0  by  passing 
through  each  point  x  =  ^  on  the  line  t  =  0  the  forward  charac- 
teristic  with  the  slope  -rr  =  u(^)  +  c(^),  and  let  u  and  c  be 

_ 

v  'We  do  not  admit  here  the  equality  sign;   the  subsequent 
results  remain  valid,  however,  also  for  a  point  of  contin- 
uity if  the  derivatives  -r—  and  -r^  become  infinite  in  an 
appropriate  manner. 
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constant  along  it.  We  set  accordingly 


(2) 


=  {u(^)  +  c(^)   t  +  ^  ,  u  =  u(^),  c  =  c(^) 


Owing  to  the  compressive  character  of  the  discontinuity 
we  have  uQ  +  cQ  <  u^  +  c-,  ,  and  hence  the  two  wave  regions  in 
the  (x,t) -plane  obtained  from  the  two  parts  x  >  0  and  x  <  0 
overlap.  The  shock  line  is  now  to  be  drawn  in  the  common 
region  in  an  appropriate  way  such  that  the  mutilated  wave 
regions  fit  together. 


d)(o) 


—  * 


Pig.  3 

Two  forward  facing  simple  waves 
separated  by  shock. 


4.  How  to  determine  the  path  of  the  shock  is  best  seen 
in  case  the  state  of  the  gas  ahead  of  the  shock  front  is  con« 
s tant :   u  =  u  =  0,  c  =  c  for  x  >  0,  initially,  and  hence 
also  later  on.  Relation  (1)  then  entails 

(3)    u  =  (1  -  ,e?)(u  +  c  -  cQ),  c  =  cQ  +  ^(u  +  c  -  cQ). 


The  relationship  between  shock  velocity  U  and  gas  velocity  u 
behind  the  shock  is  exactly  given  by 
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(4)  U  -  (1  -  y|)(U  -  C^/U)  , 

see  Manual  III>Art36(4$,  or  using  (3),  by 

/c-\  u+c   n  _  U    co 

(5)  -  1  _    - 

o        o 

Since  the  velocity  behind  the  shock  front  is  a  given  function 
u(^)  of  the  initial  parameter  ?  we  can  determine  D  from  (3) 
or  (5)  as  a  function  U(^)  of  X  .  Thus  we  have  for  the  shock 

te)  w-'W  • 

We  now  consider  the  quantity  t  along  the  shock  line  as 
function  of  £  and  express  x  as  function  of  ?  by  (2). 
Differentiation  of  equation  (2)  with  respect  to  "5  yields 

■ 

or 

(V)         Ju-u-c|£-  Ju.  +  c.(  t  ■  1  . 


{'  —  • ]S\-  {»,♦•>} 


This  is  an  ordinary  linear  differential  equation,  the  "shock 
differential  equation",  for  t  as  function  of  J   ,   which  yields 
a  unique  solution  by  virtue  of  the  initial  condition  t  =  0 
for  T  =  0.   Thus  the  propagation  of  the  shock  is  determined. 
Introducing  the  function 


(8) 


/^d(m-c) 


?=?<?>  -/'» 


with  any  convenient  constant  of  integration,  the  solution  can 
be  written  in  the  form 
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Insertion  of  this  relation  into  (2)  yields  x  as  function  of 

?  •  Thus  the  path  of  the  shock  front  is  given  in  parametric 
representation. 

For  the  numerical  treatment  various  simplifications  are 
useful  to  evaluate  the  integrals.  Chandrasekhar,  who  has 
developed  the  method  in  connection  with  a  special  case,  uses 
the  quantity  U  instead  of  F  as  parameter.  We  found  it  more 
convenient  to  employ 

u+c-c^ 

(10)  5-=  — r— 2 

C0 

as  parameter  and  to  set 

(11)  U  =  cQ  |l  +£6-  +  £    s2}    , 

a  formula  which  agrees  with  (5)  up  to  terms  of  third  order 
as  is  easily  verified  (cf.  Appendix  AI).   (Note  that  the 
dependence  of  U/c0  on  (u  +  c)/cQ  through  (5)  or  (11)  and  the 
dependence  of  £  on  U/c0  is  independent  of  y»)   Inserting  (11) 
into  (8)  we  find 

(12)  <  =  21og|-5     ,      •"*  =   (ij?)2 

whence  by   (9), 

(13>  0 1 « 8(^)2  r°  s±s^  . 

(See  Appendix  AI.) 

5.   In  the  two  problems  of  formation  and  decay  mentioned 
in  the  introduction,  we  just  have  the  situation  that  the  state 
ahead  of  the  shock  front  is  constant.  Hence  these  two  problems 
can  be  treated  by  the  method  outlined.   The  results  for  these 
two  problems  are  derived  in  Appendices  All  and  AIII.  We  give 
here  a  brief  account. 

In  the  wave  produced  by  a  pis ;on  moving  with  constant 
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1    2 

acceleration,  x  =«  bt  ,  a  shock  develops  at  the  time 

"" ™^"— """" ~~™— ~ ~"—"— ™  c  c 

ts>=  (1  -  -a)  —  »  at  the  position  x,o  =  (1  -zt)  -r~  »  (see 

Pig,  1).   The  path  of  the  shock  is  (cf.  Appendix  All),  given 
by 


r                                  3     b(t-tf)~ 
(14)   x  -  x  =  c  \{t   -  t  )  +  2-a + 


I2 

8(1-^) 


9       b2(t-t^,)3 


128  (l-^)53   "7^" 


J- 


One  observes  that  the  acceleration  of  the  shock  front  is 
initially  about  3/4  of  the  value  b  of  the  acceleration  of 
the  piston,  but  then  increases. 

6.  For  the  decaying  shock  wave  (cf •  Appendix  AIII)  we 
obtain  more  striking  results.  Suppose  the  piston  is  first 
moved  with  a  velocity  u-,  and  stopped  at  a  time  tR  at 
x  =  Xrj  =  u-,tR,  (see  Pig,  2),  A  shock  starts  moving  with 
constant  speed  at  t  =  0  from  the  piston;   it  is  overtaken  by 
a  rarefaction  wave  starting  at  the  time  tR  when  the  piston  is 
arrested.  At  the  time  t  =  t-,  of  overtaking  the  interaction 
begins  and  the  theory  of  Sec,  3  and  4  is  to  be  applied  from 
then  on.   The  motion  of  the  decaying  shock  is  found  to  be 
given  by 

(15)         t  =  tR  +  (tx  -  tR)  k2  (^)2  ,  61  >  S>   0  , 

x  =  Xp  +  (1  +  <5)  cQ  (t  -  tR)  , 

whence  asymptotically,  for  large  times  t, 


(16)  x«xR  +  c0[t  -  tQ  +  4^  y(t1-tR)(t-tR)   -  4k2(t1-tR)]  . 
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(The  values  of  the  quantities  t-,  and  k,  =  6\,/(4- S*-,)  are 
given  in  the  Appendix.)   The  width  of  the  shock  wave,  i.e. 
the  distance  between  the  shock  front  and  the  tail  of  the 
rarefaction  wave  is  then  immediately  found  to  be  asymptoti- 
cally, 

(17)   d(t)  =  co  «T(t  -  tR)  =  c0  ^4k1  ^(t^Ht^)'  -  4k^(t1-tR)J, 

Thus  it  is  seen  that  the  process  of  interaction  goes  on 
indefinitely  and  that  the  width  of  the  wave  zone  increases 
only  like  "^t-t  ' and  not  like  t  -  tR  as  the  width  of  the 
undisturbed  rarefaction  wave. 

The  shock  strength  may  be  characterized  by  the  excess 
pressure  ratio  which  is  asymptotically  given  by 


P-P0 
(18)         —2 


-  (1  +  yl)    \&    +  |  62    +  ...  j 


~(i+li2)   ^"^2  +  4^^ 


"^s"+  4k2  !i^2 

'R      1  t-tR 


Thus  we  see  that  the  shock  strength  decreases  like 

const./  T/'t-'fcR.     Tile  Pressure  distribution  over  the  wave  is 

asymptotically  given  by  the  quadratic  distribution 

(i9,         £2°~u+/f)  {e  +|  62} 


Po 


~a+P   j^^T-1+'( 


1   ^""^R       2  / 


^""'Chandrasekhar,  who  found  this  result  in  a  special  case, 
conjectured  that  this  would  be  true  in  general. 
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for  0  <  x  -  Xp  -  cQ(t  -  tR)  <  d(t),  (see  pig.  4). 
the  distribution  approaches  a  linear  one. 


Eventually 


**+$,(*-£*)   - 


Pig.  4 

Pressure  distribution  in 
decaying  shock  wave. 


7.   The  method  presented  enables  one  quite  generally 
to  handle  the  problem  of  the  flow  resulting  from  any  piston 
motion.   Of  particular  interest  is  the  case,  in  which  the 
piston  is  first  pushed  into  the  gas,  then  retracted,  and 
finally  arrested  when  it  reaches  the  original  position;  (cf. 
Appendix  AIV).   In  this  case  the  wave  zone  is  eventually 
bordered  by  two  shock  fronts.  For  this  reason  such  a  wave 
motion  has  been  termed  "N-wave".   Through  the  head  shock, 
which  was  sent  out  when  the  piston  was  first  set  into  motion, 
the  pressure  is  raised;   through  the  subsequent  rarefaction 
wave  the  pressure  drops  to  below  the  atmospheric  value,  to 
which  it  is  raised  again  through  the  tail  shock,  which  was 
sent  out  when  the  piston  was  finally  arrested. 
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Pig.  5 

Decaying  N-wave  with  head  and  tail 
shock  due  to  retraction  of  piston* 


The  width  of  the  wave  zone  in  this  case  is  asymptoti- 
cally given  by 


(20) 


d(t)    =  d1(t)    +  d2(t)    , 


with 


(21) 


V*>  =  4co  {^Nv^  t-tR1  -  ^VV} 

d2(t)  =  4c0  {-vivV^-V +  ^wj 


(note  that  kg  <  0  and  hence  dp(t)  >  0,  see  Appendix  IV); 
evidently  d(t)  is  growing  like  "yt-tR*  The  asymptotic 
pressure  distribution  is  again  given  by  (19),  but  extends 
over  the  wave  zone 


(22) 


-d^t)  <  x  -  Xj^  -  cQ(t-tR)  <  d2(t)  . 
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Htcid  Shock 


+-X 


Fig, 


Pressure  distribution  in  decaying 
N-wave  with  head-  and  tail  shock. 


The  notion  of  "N-wave11  was  introduced  recently  in  a 
paper  by  DuMond,  Cohen,  Panofsky,  and  Deeds *  ',  who  treated 
their  propagation  from  a  different  approach.  They  derive 
formulas  for  the  rates  of  change  of  the  pressure  in  the  wave 
and  at  the  shock  fronts  and  of  the  width  of  the  wave,  assuming 
that  the  wave  is  "balanced",  i.e.  that  the  product  of  the 
pressure  ratios  of  head  and  tail  shock  equals  one.   (See  (6), 
(7),  (8)  in  their  paper;   in  (6)  and  (7)  a  factor  Y  in  the 
denominator  was  omitted.)   These  formulas  are  in  agreement 
with  the  first  order  contributions  to  our  formulas  (19),  (20). 
(21),  from  which  they  could  be  derived. 


vn 


See,  "A  Determination  of  the  Wave  Forms  and  Laws  of  Propa- 
gation and  Dissipation  of  Ballistic  Shock  Waves"  by  J.  W.  M. 
DuMond,  E.  R.  Cohen,  W.  K.  H.  Panofsky,  and  E.  Deeds,  re- 
printed from  The  Journal  of  the  Acoustical  Society  of  America, 
Vol.18,  No.  1,  97-118,  July,  1946.   See  also  L.  D.  Landau, 
J.  Phys.  Acad.  Sci.  U.S.S.R.,  Vol.  6,  229-250,  1942.   In  a 
paper  by  B.  Cassen  and  J.  Stanton,  Theory  of  Non-Steady  Shock 
"waves,  of  which  an  abstract  appeared  in  the  Bulletin  of  the 
American  Physical  Society  V«ol.  22,  No.  2,  May  1947,  similar 
problems  are  discussed  by  an  approach  apparently  more  closely 
related  to  ours. 
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8.  One  naturally  asks  how  much  faith  may  be  put  In  the 
asymptotic  representations  obtained.   A  few  critical  remarks 
are  called  for.   The  method  employed  can  be  considered , as 
giving  the  correct  terms  up  to  the  second  order  of  an  expan- 
sion with  respect  to  powers  of  the  shock  strength.   Therefore 
the  position  of  the  shock  front  at  a  definite  time  t  is  given 
with  arbitrarily  good  accuracy  if  the  strength  of  the  shock 
initially  is  made  sufficiently  small.   One  may  be  tempted  to 
infer  from  this  fact  that  our  formulas  for  the  position  x  of 

a  shock  front  in  a  given  problem  represent  the  actual  position 
with  increasing  accuracy  for  increasing  time  t,  since  the 
shock  strength  approaches  zero  as  the  time  increases.   This 
inference  would,  however,  not  be  correct.   For  the  case  of  the 
decaying  shock  wave  treated  in  Sec.  6  this  would  mean,  for 
example,  that  the  coefficient  k,  in  formula  (16)  was  indepen- 
dent of  the  approximation  method,  which  cannot  be  expected  to 
be  true.   All  that  can  be  expected  is  that  an  asymptotic 
representation  of  the  type  (16)  holds  with  coefficients  not 
differing  much  from  those  in  formula  (16).   Thus  the  deviation 
of  the  actual  position  x  of  the  shock  from  the  one  given  by  (16) 
must  be  expected  to  increase  like  a  quantity  <*  JTf   but  on  the 
other  hand  the  factor  «*•  may  be  expected  to  be  small. 

In  order  to  determine  the  accuracy  of  the  method  one  may 
refine  the  description  of  the  shock  wave  by  including  terms  of 
third  order  in  the  shock  strength.   Results  of  such  an  investi- 
gation will  be  reported  in  a  later  issue. 

9.  For  the  treatment  of  the  problem  of  a  shock  front  separ- 
ating two  simple  waves  (neither  of  which  is  a  constant  state)  (cf. 
Fig.  3)  the  method  explained  in  Sec.  4  is  no  longer  applicable. 
According  to  the  argument  of  Sec.  4,  we  may  again  assume  that 

up  to  second  order  the  simple  waves  are  not  influenced  by  the 
shock;  but  it  does  not  seem  possible  to  determine  the  motion 
of  the  shock  front  from  a  linear  differential  equation.   Of 
particular  interest  is  the  case  in  which  the  shock  begins  with 
initial  strength  zero;   such  a  s:  .uation  arises,  for  example, 
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if  a  shock  develops  within  a  simple  wave.   In  that  case  an 
expansion  with  respect  to  powers  of  the  time  will  give 
sufficient  information  about  the  early  stages  of  the  develop- 
ment of  the  shock,  (cf.  Appendix  AV). 

Let  us  assume  the  simple  wave  to  be  given  by  one  repre- 
sentation 

(23)  x  =  »  t  +  A(«  )  , 
with 

(24)  °  =  u  +  c  ; 
u  and  c  are  then  given  by 

(25)  u  =  u,  +  (1  -  /£)(*>  -  af),    c  =  c,  +   ^(«-  «,}, 

u»  +  %       =  a9 

with  appropriate  constants  u^  and  c,  •  Let  us  assume  that 
the  function  A(«  )  describing  the  initial  distribution  of u 
satisfies  the  condition 

'(26)  A«(«»  )  <  0  , 

except  for  Cj    =    &p     where 

(27)  A'(^)=A"(^)=0,  A'"(^)<0   . 

Then  the  straight  forward  characteristics  form  an  envelope 
given  through  the  parameter  representation 

(28)  t=-A'(w),x=A(0)-0A'(0) 

with  a  cusp  at  t  =  0,  x  =  0.  Consequently,  a  shock  will 
develop.  As  will  be  shown  in  the  Appendix  AV,  the  expansion 
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of  tiie  representation  of  this  shock  begins  with  the  terms 
(29)        x  =  C^t  +  ^(1  +  k"1^  )c^t2A  > 
in  which 

(so)    k  =  -  \  c2  a"'(  <y,),  ^4C^3  A'n,(  %>  • 

It  is  remarkable  that  up  to  the  order  considered  the  shock 
depends  only  on  the  initial  distribution  of  CJ    =  u  +  c  and 
on  the  local  value  c  ^  of  c  but  not  explicitly  on  the  value 
of  y. 

These  formulas  can  be  used  to  determine  the  formation 
of  a  shock  in  a  compression  wave  which  is  produced  when  a 
piston  is  pushed  into  the  gas  beginning  with  acceleration 
zero.   This  will  be  explained  in  detail  in  Appendix  AVI. 

10.   The  decay  of  steady  two-dimensional  shock  waves 
can  be  treated  in  much  the  same  way  as  the  decay  of  moving 
one-dimensional  shock  waves  has  been  treated.  The  reason 
is  that  for  steady  two-dimensional  waves,  just  as  for  one- 
dimensional  unsteady  waves,  the  transition  relations  for 
shocks  and  simple  waves  are  the  same  up  to  terms  of  second 
order. 

Consider  a  backward-facing  simple  wave  that  is  forward 

inclined.  Suppose  further  that  the  incoming  supersonic  flow 

with  the  constant  velocity  u=q  >  c  ,  v  =  0  is  turned  into 

J  ^o    o 

a  flow  with  velocity  u  =  q  cos  ©,  v  =  q  sin  0,  (see  Pig.  7). 
The  strength  of  this  wave  may  be  characterized  by  the  angle 
©  and  the  changes  of  all  quantities  across  the  wave  may  be 
expanded  in  powers  of  0.  Similarly,  consider  a  backward- 
facing  forward-Inclined  shock  front,  (see  Pig.  8)  across 
which  the  flow  with  velocity  (q  ,0)  is  turned  into  a  flow 
with  velocity  u  =  q  cos  9,  v  =  q  sin  ©.   Then  again  the  changes 
of  all  quantities  can  be  expanded  in  powers  of  ©.   This  ex- 
pansion now  agrees  with  the  expansion  for  the  simple  wave  up 
to  terms  of  second  order,  (see  Appendix  BII).  For  this  reason 
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we  may  in  second  approximation  consider  the  flow  behind  the 
shock  as  if  it  had  resulted  through  a  simple  wave. 


Pi 


<=>• 


Backward-facing  forward- 
inclined  simple  wave  in- 
dicated through  a  set  of 
straight  Mach  lines. 


Pig.  8 


Backward-facing  forward- 
inclined  shock.  Forward- 
inclined  Mach  lines  are 
shown  on  both  sides. 


This  fact  may  also  be  expressed  as  follows:   consider 
in  the  (u,v) -plane  the  shock  polar  and  the  epicycloid  through 
the  same  point  (q  ,0)  employing  the  same  value  for  the  critical 
speed  cv,,  (see  Manual  III,  35,  and  IV,  61,  69)  „     Then  these 
two  curves  agree  up  to  second  order  in  the  sense  that  their 
representations  in  terms  of  the  parameter  ©  agree  up  to  second 
order.   If,  therefore,  the  flow  velocity  u,v  behind  a  shock 
with  a  certain  angle  0  of  flow  direction,  represented  by  a 
point  on  the  shock  polar,  is  replaced  by  a  velocity  with  the 
same  angle  0  represented  by  a  point  on  the  epicycloid,  then 
the  shock  relations  for  u  and  v  are  satisfied  up  to  terms  of 

second  order  in  the  shock  strength.   The  sound  speed  c  is 

2  2         2   2    2 

determined  from  Bernoulli's  law  juQ.     +    (1  -  /*)c  =  c„  and 

p  and  p  are  functions  of  c,  if  we  assume  the  entropy  to  be  the 
same  on  both  sides. 

Suppose  now  a  supersonic  -low  along  a  straight  wall  with 
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a  bump  is  to  be  determined,  (see  Fig.  9).  We  may,  for  ex- 
ample, consider  the  upper  part  of  an  airfoil  cross  section 
as  such  a  bump.  We  assume  that  the  incoming  flow  is  super- 


Pig,  9 

Plow  along  a  straight  wall  with  a 
bump.  The  resulting  shock  and  for- 
ward-facing Mach  lines  are  shown. 


sonic  M  =  Q0/c  >   1>  that  the  profile  of  the  bump  begins 
with  a  sharp  angle  9,  of  inclination  against  the  straight 
wall,  which  we  take  as  x-axis,  and  that  9-,  is  so  small  that 
the  turning  of  the  flow  through  the  angle  9-,  may  be  achieved 
by  a  shock  starting  at  the  tip*  The  continuation  of  this 
shock  varies  in  strength  and  inclination  (unless  the  bump 
is  straight);  hence  the  flow  behind  the  shock  front  is  ro- 
tational and  carries  a  non-constant  entropy.   In  second  order, 
however,  this  flow  is  just  the  simple  wave  flow  that  would 
have  resulted  if  the  flow  had  turned  through  the  angle  9,  by 
means  of  a  simple  compression  wave.  This  simple  wave  flow 
corresponds  in  the  (u,v) -plane  to  the  epicycloid  which  begins 
at  the  point  (q  ,0)  and  carries  the  same  entropy  as  the  flow 
ahead  of  the  shock.   Once  this  simple  wave  has  been  determined 
we  shall  obtain  approximately  the  shock  line  as  the  solution 
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of  an  appropriate  linear  differential  equation.   This 
equation  is  based  on  an  expression,  valid  up  to  terms  of 
second  order,  for  the  slope  of  the  shock  line  in  terms  of 
the  states  on  both  sides.   This  shock  line  cuts  off  a  part 
of  the  full  simple  wave  flow.   The  flow  so  described  satis- 
fies all  conditions  correctly  up  to  the  second  order  since 
the  states  on  both  sides  of  the  shock  line  are  represented 
in  the  (u,v) -plane  by  points  on  the  same  epicycloid  and  have 
the  same  entropy. 

In  order  to  determine  the  simple  wave  flow  that  repre- 
sents the  flow  behind  the  shock  line  in  second  approximation, 
we  suppose  that  the  bump  profile  is  given  by  a  representation 

(31)     y  =  Y(x),  x  >  0;   Y(0)  =  0   . 

Through  every  point  x  =  f  of  the  profile  a  straight  Mach 
line  passes,  given  by  the  equation 

(3^)  y  =  t(x  -  £  )  +  Y($) 

in  which  the  slope  t  of  the  Mach  line  is  a  given  function 
t  =  t(9)  of  the  angle  of  inclination  0  of  the  flow  with  the 
x-axis.  Since  the  gas  moves  along  the  profile,  this  angle  9 
is  connected  with  the  function  Y(£ )  through 

(33)  tan  ©  =  Y'($)  . 

The  slope  t  is  given  by  the  relation 

(34)  t  =  tanK  +  9)  , 

in  which  oC  is  the  Mach  angle,  connected  with  the  Mach  number 
M  =  q/c  through 


(35) 


M  =  l/sinof  or  1/  7M2-1      =  tan  «<   . 
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The  Mach  number  M  in  turn  can  be  expressed  through  the  flow 
speed  q  and  the  connection  between  the  flow  speed  q  and  the 
angle  0  of  flow  direction  can  be  read  off  from  the  epicycloid 
through  the  point  (q  ,v),   In  this  way  t  can  be  determined 
as  function  of  9. 

To  determine  the  shock  line  we  may  employ  any  relation 
which  expresses  the  slope 

(36)  s  =  tan  p 

of  the  shock,  or  its  angle  of  inclination  p,  in  terms  of 
quantities  on  both  sides  of  it.   It  is  natural  to  use  a 
relation  connecting  s  with  the  angle  ©.  Along  the  shock  line 
we  have 


(37)  g 


ay  _ 


s 


• 


Inserting  s  as  function  of  0,  and  ©  as  function  of  ^  ,  into 
the  relation 

derived  from  (32),  we  obtain  the  linear  differential  equation 

(38)  (t  -  s)§|  +  £|x  =  t  +f  £|  -  Y»(J)  , 

for  x  as  function  of  I      ,  which  is  to  be  solved  under  the 
initial  condition 

(39)  x  =  0  for   ^   =  0  . 

For  the  solution  of  this  differential  equation  it  is 
convenient  to  expand  t  in  powers  of  the  angle  ©  and  to  retain 
only  the  terms  up  to  the  second  order.  We  then  obtain,  (see 
Appendix  Bill) 


with 


(40)         t  =  tQ  +  a-^  +  a2to1©2  +  ... 


(41)  a,  .JLyU  +  tJ)8 


We  note  that  t  is  connected  with  the  Mach  number  M  of  the 
incoming  stream  through 

(42)    M;2  =  t2(l  +  t^)"1  or  t2  =  M^2(l  -  m;2)"1   . 


We  further  expand  s  in  powers  of  t  -  t  ,  obtaining 

(43)       s  =  t0  +  |(t  -  tQ)  ^^(t  -  t0)2 

with 

l-3t 

(44)         k='I^£=^-M°2   • 

O 

Equation  (38)  then  takes  the  form 

(t-t0)-lrt^(t-t0^ 


^-— T2^  <**-«?» 


(t-t^-kt^ft-t,,) 


The  solution  of  this  equation  is 


(46) 


x(£)    =  2(- 
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l-kt"1(t-t^) 


t-i 


■)' 


/  ^ *   d[  Jf  't'-Y(l')] 

Jo  [l-kt^tt* -t  )]3 


in  which  t  =  t(®  (£))  and  t«  =  t(  ®  (£»));   y(g)  is  then 
obtained  by  insertion  in  (32). 

A  general  conclusion  can  be  drawn  from  this  formula. 
When  the  parameter  £  approaches  a  value  £  =  *m  for  which 
0  approaches  zero  and  hence  t  approaches  t  the  factor  in 
front  of  the  integral  becomes  <d  ,  and  thus  both  x  and  y 
become  infinite.   In  other  words  the  shock  front  penetrates 
only  into  that  part  of  the  simple  wave  which  Issues  from  the 
section  of  the  bump  between  its  beginning  £  =  0,  and  its 
peak,  £  =  x   .   The  shock  line  extends  to  infinity  without 
ever  reaching  the  Mach  line  which  Issues  from  this  peak  and 
asymptotically  it  becomes  parallel  to  this  Kach  line. 


X.m 


Pig.  10 

The  shock  front  will  asymptotically 
become  parallel  to  the  Kach  line 
issuing  from  the  psak. 
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To  investigate  more  closely  the  asymptotic  behavior  of 
the  shock  front  at  great  distances  from  the  axis  we  expand 
all  quantities  with  respect  to  negative  powers  of  y.   we 
introduce  the  quantity  $  **'  by 

(47)   ef  =  2t;2  f  (t-t0)[l-l:t;1(t-t0)]-3d(^  t-Y(£)) 

^o 

and  find  from  (46) 

(48)  x  =  cr{t2(t-tcr2  -  2kt0(t-tor1  + ...  , 

from  which 

(49)  t-tQ  =  tQ  («T  VV1^  .  k^x-l  +  ###  m 
Inserting  relation  (49)  in  (32)  gives 

(50)  y  =  ym  +  ^(xrfl^)  +  tQ  ^x1/2  -  tQk  g  +  . . . 
Expressing  x  in  terms  of  y  we  find 

(51)  x  =  x^  +  t^y-yj  -  d(y)  , 
in  which 

(52)  d(y)  =  6 !/2t-l/2yl/2  _  (1  +  k)  ^  +  ###  % 


'"""'If  the  profile  is  given  by 

y  =  7»(2 1 "  (t,2) 

we  have, 

S  "  ^l^m  +  ^ar2al+6kal+2a2^o2ym  +   ••'      ' 
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This  quantity  d  =  d(y)  is  the  distance  the  shock  front  stands 
ahead  of  the  Mach  line  x  =  x  +  t"  (y-y) issuing  from  the  top 
of  the  profile;  we  see  that  this  distance  increases  like  the 
square  root  of  the  distance  y  from  the  axis . 

The  strength  of  the  shock  is  in  first  order  proportional 
to  t-t  ;   hence  the  strength  of  the  shock  decreases  like  the 
reciprocal  of  the  square  root  of  the  distance  from  the  y-axis. 
To  calculate  this  strength  as  a  function  of  y  one  may  proceed 
as  follows.   One  inserts  (51)  into  (49),  and  (49)  into  (43); 
then  one  uses  the  formula 

2      2 

P-P  2        2  n+t 

(53)  Ho=Ji±S-±o      § 

p   q2Q  1+S2  1+tJ 

(cf.  Manual  (52)  IV-69,  p.  193).   The  expansion  of  the  rela- 
tive pressure  increase  (p-p  )/p  with  respect  to  negative 
powers  of  y  '  *'  could  easily  be  derived.  The  term  of  first 
order  obtained  this  way  is 

t*±\  P"P°   n,  2,  to  &   v  -1/2 

(54)  —  =  (1-f  h)       — g   y     +  ... 

o 


Also  the  gas  velocity  behind  the  shock  and  the  flow  angle  9 
are  easily  expressed  in  terms  of  s  =  tan  p,  (cf.  (50),  (52), 
IV-69,  p.  192,  193).   The  angle  0  as  a  function  of  y  could  be 
more  directly  obtained  by  inverting  relation  (40). 

The  critical  remarks  made  in  Sec.  8  apply  in  an  obvious 
manner  to  the  problem  treated  in  this  article.  Dumond,  Cohen, 
Panofsky  and  Deeds  (see  p.  13  )  have  made  the  very  interesting 
remark  that  the  width  of  the  shock  wave  is  determined  asymptoti- 
cally solely  by  the  shock  strength  and  the  distance  from  the 
profile,  independently  of  the  specific  shape  of  the  profile. 
We  find  the  same  result  from  (52)  and  (54)  in  first  order: 
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P^£  =  (1  +   2)  J°  Oil 


(55)  y!'  (1  ♦  /f)  7^3 

'o 


The  coefficient  here  differs,  however,  from  the  coefficient 
given  in  that  paper  by  the  factor  2tQ/  v-L+t   =  2  sin  «C  (see 
I.e.  (9)). 
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